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Abstract 

Computational  fluid  dynamics  analysis  was  employed  to  investigate  the  performance  of  proton  exchange  membrane  fuel  cells  (PEMFCs)  with 
different  channel  geometries  at  high  operating  current  densities.  A  3D,  non-isothermal  model  was  used  with  a  single  straight  channel  geometry. 
Both  anode  and  cathode  humidifications  were  included  in  the  model.  In  addition,  phase  transportation  was  included  in  the  model  to  obtain  the 
total  water  management  for  systems  operating  at  different  current  densities.  The  simulation  results  showed  that  a  rectangular  channel  cross-section 
gave  higher  cell  voltages  compared  with  trapezoidal  and  parallelogram  channel  cross-sections.  However,  the  trapezoidal  channel  cross-section 
facilitated  reactant  diffusion,  leading  to  more  uniform  reactant  and  local  current  density  distributions  over  the  reacting  area,  and  thus  to  a  lower 
cathode  overpotential  of  the  cell.  Simulations  of  the  three  different  channel  cross-sections  using  the  same  boundary  conditions  showed  that  among  the 
cell  geometrical  parameters,  the  shoulder  width  is  one  of  the  most  influential  in  terms  of  its  impact  on  cell  performance.  Simulations  using  different 
channel-shoulder  width  ratios  showed  that  at  high  operating  current  densities,  Ohmic  losses  significantly  increase  with  decreasing  shoulder  width. 
In  contrast,  a  smaller  shoulder  width  facilitates  the  distribution  of  reactants  and  helps  to  reduce  concentration  losses.  The  simulations  disclosed  the 
existence  of  an  optimum  channel-shoulder  width  ratio  that  gives  the  highest  cell  voltage  under  high  current  density  operating  conditions.  Under 
such  conditions,  however,  the  cell  performance  deteriorated  dramatically  with  decreasing  shoulder  width,  even  when  higher  reactants  flow  rates 
and  inlet  velocities  were  used. 
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1.  Introduction 

The  proton  exchange  membrane  fuel  cell  (PEMFC)  is  a 
promising  technology  for  various  applications,  including  auto¬ 
mobiles  and  stationary  power  plants.  However,  PEMFCs  require 
further  improvements  if  they  are  to  be  used  commercially  in  the 
automobile  industry,  particularly  in  regard  to  their  power  den¬ 
sity  [1,2].  To  obtain  high  power  density,  voltage  losses  of  the  cell 
must  be  minimized  and  the  cell  must  be  operated  at  high  operat¬ 
ing  current  density.  It  is  well  established  that  the  major  source  of 
voltage  loss  in  PEMFCs  is  the  oxygen  reduction  reaction  (ORR) 
at  the  cathode  side  [3,4].  Furthermore,  the  water  produced  at  the 
cathode  side  leads  to  concentration  losses  under  high  current 
density  operation.  Many  researchers  have  highlighted  this  issue 
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and  emphasized  the  need  for  optimization  of  geometrical  and 
physical  parameters  along  with  proper  water  and  thermal  man¬ 
agement.  Results  obtained  using  computational  fluid  dynamic 
(CFD)  analysis  have  indicated  the  need  to  reduce  concentration 
losses  in  order  to  improve  cell  performance  [4-10]. 

Hydrogen  and  oxygen  diffusion  through  the  gas  diffusion 
layer  (GDF)  and  the  distributions  of  these  species  at  the  reacting 
area  are  also  important  determinants  of  fuel  cell  performance. 
One  of  the  easiest  ways  to  facilitate  the  flow  of  reactants  to 
the  reacting  area  is  to  change  the  geometrical  configuration  of 
the  channel  (e.g.,  length,  width,  height  and  shoulder  (bipolar 
plate)  width)  and  the  channel  flow  pattern.  A  number  of  studies 
have  endeavored  to  improve  fuel  cell  performance  by  varying  the 
channel  configuration  [8,9,11-17].  It  has  been  suggested  that  use 
of  a  trapezoidal  channel  cross-section  rather  than  a  rectangular 
or  square  cross-section  may  improve  the  cell  performance  [18]. 
Various  serpentine  flow  field  designs  for  both  the  anode  and  cath¬ 
ode  channels  have  also  been  examined  to  obtain  higher  output 
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Nomenclature 

Acv 

specific  surface  area  of  the  control  volume  (m-  1 ) 

Areacv 

surface  area  of  the  control  volume  (m2) 

C\va 

concentration  of  water  vapor  at  the  anode 
(mol  m-3) 

CWc 

concentration  of  water  vapor  at  the  cathode 
(mol  m-3) 

DH2,1 

diffusion  coefficient  of  H2  in  a  liquid  water  film 
(6.3  x  10-9  m2  s-1) 

£>02,i 

diffusion  coefficient  of  O2  in  a  liquid  water  film 
(2.4  x  10-9  m2  s-1) 

£>w 

diffusion  coefficient  of  water  (m2  s-1) 

F 

Faraday  constant  (96,487  C  mol-1) 

4xn 

enthalpy  of  water  formation  (kJ  Kmol-1) 

hfo 

enthalpy  of  vaporization  of  water  (kJkg-1) 

hh2.\ 

Henry’s  law  constant  for  H2  in  a  liquid  water  film 
(8.9  X  109  Pa) 

H(h.\ 

Henry’s  law  constant  for  O2  in  a  liquid  water  film 
(2.12  x  1010  Pa) 

I 

local  current  density  (Am-2) 

4vg 

average  current  density  (Am-2) 

exchange  current  density  for  reaction  K  (A  m-2) 

L 

length  (m) 

37m.  dry 

equivalent  weight  of  a  dry  membrane  (kg  mol-1) 

Mn 

molecular  weight  of  species  n  (kg  mol-1) 

%,1 

mass  fraction  of  liquid  water 

mass  n 

mass  of  species  n  (kg) 

n 

no.  of  electrons 

nd 

electro-osmosis  drag  coefficient 

psat 

Uv  ,k 

saturation  vapor  pressure  of  water  in  stream 
k  (Pa) 

P 

pressure  (Pa) 

Pn 

partial  pressure  of  species  n  (Pa) 

R 

universal  gas  constant  (8.314  J  (mol  K)-1) 

r 

condensation  rate  (s-1) 

S 

source  term 

k 

liquid  water  film  thickness  (m) 

hn 

membrane  thickness  (m) 

T 

temperature  (K) 

V 

velocity  (ms-1) 

Tcell 

cell  voltage  (V) 

Voc 

open  circuit  voltage  (V) 

W 

width  (m) 

xik 

mole  fraction  of  species  i  in  stream  k 

Greek  letters 

a 

net  water  flux  per  proton 

p 

permeability  (m2) 

£ 

porosity  of  the  gas  diffusion  layer 

T] 

overpotential  (V) 

X 

water  content  in  the  membrane 

M 

dynamic  viscosity  (kg  s  m-2) 

Pm,  dry 

density  of  a  dry  membrane  (kg  m-3) 

P 

density  of  the  mixture  (kgm-3) 

crm  membrane  conductivity  (S  m  1 ) 

§  stoichiometric  rate 

Subscripts  and  superscripts 
a  anode 

c  cathode 

ch  channel 

cv  control  volume 

e  electrochemical  reaction 

glif  gas  liquid  interface 

H2  hydrogen 

k  anode  or  cathode 

1  liquid 

O2  oxygen 

reac  reacting 

sat  saturated 

v  vapor 

§  dummy  variable  for  direction  x,  y  or  z 


from  PEMFCs  [19-27].  All  researchers  have  focused  on  achiev¬ 
ing  a  uniform  reaction  over  the  reacting  area  by  facilitating  the 
distribution  of  oxygen  in  the  reacting  area.  Under  high  operating 
current  density  conditions,  however,  oxygen  molecules  struggle 
to  reach  the  reacting  area,  especially  in  the  downstream  portion 
of  the  channel. 

In  PEMFCs,  the  cell  performance  also  depends  on  various 
physical  and  electrochemical  parameters.  Recently,  Tao  et  al. 
[28]  and  Min  et  al.  [29]  examined  the  parameter  sensitivity  (11 
parameters)  of  a  PEMFC  in  a  3D,  non-isothermal  model.  In  these 
simulations,  the  parameters  were  linked  so  as  to  reflect  the  situa¬ 
tion  in  real  PEMFCs;  this  meant  that  a  variation  of  one  parameter 
that  had  a  positive  influence  on  cell  performance  could  lead  to 
variations  in  other  parameters  that  negatively  affected  the  per¬ 
formance  [29] .  These  studies  showed  that  PEMFC  performance 
was  sensitive  to  variations  in  the  transfer  coefficient,  exchange 
current  density,  water  vapor  diffusion  coefficient,  oxygen  dif¬ 
fusion  coefficient,  GDL  porosity  and  membrane  conductivity. 
However,  other  parameters,  such  as  the  stoichiometric  rate,  pres¬ 
sure  and  temperature,  also  influenced  the  fuel  cell  performance 
[8,15,30,31],  especially  the  diffusion  of  reactants  to  the  reacting 
area.  From  a  geometrical  standpoint,  the  shoulder  (bipolar  plate) 
and  channel  configurations  also  influence  the  flow  of  reactants 
and  have  a  marked  effect  on  fuel  cell  performance.  Moreover, 
these  geometrical  parameters  also  affect  the  inlet  velocity  at  the 
anode  and  cathode  [32-34].  Usually  the  inlet  velocities  of  the 
reactants  into  the  channels  are  defined  by  the  following  equation: 


£k  4vg  R  2k  Areac 

nFPXi  Ach 


where  §,  /avg,  R ,  T,  n ,  F  and  P  are  the  stoichiometric  rate,  aver¬ 
age  current  density,  molar  gas  constant,  temperature,  number  of 
electron,  Faraday’s  constant  and  pressure,  respectively.  2Q  is  the 
mole  fraction  of  species  2,  Areac  and  Ac h  are  the  reacting  area  and 
channel  cross-sectional  area,  respectively.  The  subscript  k  refers 
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to  anode  or  cathode.  A  higher  inlet  velocity  facilitates  passage 
of  the  reactants  to  the  reacting  area  and  leads  to  a  more  uniform 
reaction  in  the  cell.  Many  researchers  have  studied  the  effects  of 
various  channel-shoulder  width  dimensions  on  PEMFC  perfor¬ 
mance  [8,11-14,16,27].  When  modifying  the  PEMFC  configu¬ 
rations  like  channel  width,  shoulder  width  or  reacting  area  lead 
to  change  the  boundary  conditions  like  reactants’  inlet  veloci¬ 
ties  or  mass  flow  rates.  However,  most  researchers  have  either 
not  reported  the  inlet  velocities  and  flow  rates  of  the  respective 
channels  or  have  not  adequately  defined  the  boundary  condi¬ 
tions  for  each  cases.  In  a  study  of  channel  size  effects  using 
a  2D  isothermal  model,  in  which  the  channel-shoulder  width 
ratio  was  maintained  constant  (while  varying  both  the  chan¬ 
nel  and  shoulder  size),  Guvelioglu  and  Stenger  [13]  found  that 
cell  performance  was  superior  in  systems  with  smaller  chan¬ 
nel  and  shoulder  sizes.  Natarajan  and  Nguyen  [8]  and  Lin  and 
Nguyen  [12]  used  a  2D  model  to  study  the  effect  on  cell  perfor¬ 
mance  of  a  variety  of  parameters;  specifically,  they  examined 
the  effects  on  a  cell  with  a  given  membrane  width  of  vary¬ 
ing  the  number  of  channels  and  shoulders  while  maintaining 
a  constant  channel-shoulder  width  ratio,  and  also  of  varying 
the  channel-shoulder  width  ratio.  They  found  that  oxygen  and 
liquid  water  were  transported  over  the  shoulder  area  more  eas¬ 
ily  in  configurations  with  smaller  shoulder  widths,  leading  to 
enhanced  cell  performance.  Lum  and  McGuirk  [14]  and  Sun  et 
al.  [16]  studied  the  effect  of  shoulder  width  on  fuel  cell  per¬ 
formance  in  2D  and  3D  models,  respectively,  by  varying  the 
shoulder  width  while  maintaining  a  constant  channel  width. 
They  showed  that  for  systems  with  smaller  shoulder  widths, 
more  oxygen  diffuses  to  the  catalyst  layer  through  the  GDL 
because  of  the  shorter  distance  and  that  water  transportation  is 
improved.  Moreover,  their  results  indicated  that  the  local  cur¬ 
rent  density  distribution  is  significantly  improved  by  decreasing 
the  shoulder  width.  Recently,  Chiang  and  Chu  [17]  conducted 
numerical  simulations  of  a  3D  isothermal  model  of  a  PEMFC 
for  various  channel  configurations  (channel  length,  height  and 
shoulder  width),  with  particular  emphasis  on  the  effect  of  chan¬ 
nel  height  on  cell  performance.  They  carried  out  simulations 
with  different  combinations  of  channel  width  and  height  while 
maintaining  the  same  channel  cross-sectional  area  and  constant 
flow  rate.  Their  results  indicated  that  flat  channels  (i.e.,  smaller 
height)  gave  better  performance. 

In  the  present  study,  we  performed  simulations  to  investigate 
the  performance  of  PEMFCs  at  high  operating  current  den¬ 
sity  for  various  channel  cross-sectional  configurations,  while 
maintaining  the  same  reactant  flow  rates  and  inlet  boundary  con¬ 
ditions.  In  addition,  we  examined  the  different  types  of  voltage 
losses  of  the  cell  with  related  water  and  thermal  management  at 
high  operating  current  density.  Emphasis  was  placed  on  exam¬ 
ining  the  effect  of  varying  the  channel-shoulder  width  ratio 
for  a  given  channel  height  under  high  current  density  oper¬ 
ation.  Particular  attention  was  given  to  investigating  the  cell 
performance  with  the  model  boundary  conditions  governed  by 
Eq.  (1)  and  also  to  maintain  the  same  boundary  conditions  for 
different  channel-shoulder  width  ratios.  For  each  geometrical 
configuration,  the  cell  performance  was  investigated  in  detail 
by  examining  the  reactants,  local  current  density,  membrane 


conductivity,  anode  and  cathode  overpotentials,  and  water  and 
thermal  management.  A  3D  model  with  a  single  straight  channel 
geometry  was  used  for  the  computational  analysis.  The  simu¬ 
lation  results  were  validated  by  comparison  with  results  in  the 
literature. 


2.  Numerical  simulation 


Schematic  diagrams  of  the  computational  domain  and  its 
cross-sectional  view  are  shown  in  Fig.  1.  The  system  con¬ 
sists  of  two  channels  (for  hydrogen  and  air)  between  which  the 
membrane  electrode  assembly  (MEA)  is  sandwiched.  There  are 
bipolar  plates  (also  known  as  shoulders)  on  both  the  anode  and 
cathode  sides,  which  function  as  current  collectors  with  high 
electronic  conductivity.  Humidified  hydrogen  and  air  are  intro¬ 
duced  into  the  respective  channels  in  quantities  determined  by 
the  stoichiometric  rate  and  average  current  density.  The  follow¬ 
ing  assumptions  are  made  in  the  model:  ideal  gas  mixture,  steady 
state,  laminar  flow,  isotropic  and  homogeneous  porous  gas  diffu¬ 
sion  layer,  homogeneous  two  phase  flow  and  negligible  potential 
drop  in  the  bipolar  plates. 

The  governing  equations  for  the  numerical  simulation  are 
conservation  of  mass,  momentum  transport,  species  transport 
and  energy  equations. 

Conservation  of  mass  equation: 

V  •  ( pu )  =  Sm&  +  Smc  (2) 

The  source  terms  are 


3Vn:i  —  ^H?  T  ^wvn  T  ^wln  T  $•, 


awve  > 


(2.1) 


Nnc  —  ^02  3“  ^wvp  T  Swlp  T  ^cwve  j 


(2.2) 


where 


4Th9  Acv/ 


Ac\I 


SH?  = - l  ^  and  So2  = - 2A^'  .  (2.3) 


IF 


4  F 


The  mass  contributions  of  the  water  vapor  at  the  anode  and 
cathode  sides  are 


^awvc  — 


^h2o  ACYaI 
F 


(2.4) 


^CWVe  —  ^CWVt  +  ^CWVp  — 


47r9o Acv/  .  cxMh7oAcv/ 


IF 


+ 


F 


Fig.  1.  Single  straight  channel  flow  field  for  a  PEMFC.  (a)  Geometry  and  (b) 
cross-sectional  view. 
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(1  +  2a)M}i2oAcvI 
IF 


(2.5) 


The  phase  change  between  water  vapor  and  liquid  water 
depends  on  the  partial  pressure,  and  is  defined  as 


ofvmasS^  of v/ of  v ) 

(1  -  (PS*/P)) 


X 


x  r. 


(2.6) 


r  is  the  condensation  rate  and  set  as  1  s_1  in  the  model.  Eq.  (2.6) 
shows  the  way  of  phase  change  in  water  by  checking  the  partial 
pressure  of  water  vapor  compared  with  saturated  pressure.  If  the 
partial  pressure  of  water  vapor  is  greater  than  the  saturated  pres¬ 
sure,  this  water  vapor  will  condense  to  liquid  water.  If  the  partial 
pressure  of  water  vapor  is  lower  than  the  saturated  pressure,  the 
liquid  water  will  evaporate  to  water  vapor. 

In  Eqs.  (2.4)  and  (2.5),  a  is  the  net  water  transfer  coefficient 
per  proton  through  the  membrane  and  defined  as 

FDw[Cwc  —  Cwa] 

a  =  nd - - -  (2.7) 

ttm 

The  electro-osmotic  drag  coefficient  (^d)  and  water  diffusion 
coefficient  (£>w)  can  be  correlated  with  membrane  water  content 
(A)  [35]  as  follows: 

«d  =  0.0029A2  +  0.05A  -  3.4  x  1(T19.  (2.8) 


Momentum  transport  equation: 

V  •  (pww)  —  — V P  T  V  •  (/x  Vw)  Spj .  (3) 


where  Spj  is  the  sink  source  term  for  porous  media  in  the  x,  y  and 
z  directions.  As  the  pressure  drops  in  porous  media,  Darcy’s  law 
has  taken  into  account  in  the  model.  The  source  term  is  defined 
as 


(3.1) 


Here,  /3  is  the  permeability. 

General  species  transport  equation: 


V  •  ( pmnu )  =  V  •  ( Jn )  +  Ss. 


where  n  denotes  H2,  O2,  water  vapor  or  liquid  water.  The  source 
terms  are  the  same  as  those  for  the  conservation  of  mass  equa¬ 
tion.  The  transportation  of  the  species  is  solved  with  bulk  mixture 
velocities  and  with  diffusion  mass  fluxes.  The  diffusion  mass 
flux  of  each  species  is  evaluated  with  binary  diffusion  coef¬ 
ficients  [36]  and  is  reduced  by  50%  in  the  diffusion  layer  to 
account  for  the  effect  of  porosity  and  tortuosity  of  the  pores 
[37].  It  is  assumed  that  the  liquid  water  is  in  small  droplet  forms 
and  freely  suspended  in  the  gas  flow.  This  allows  that  two  phases 
flow  together  (homogeneous  two  phase  flow)  at  normal  operat¬ 
ing  conditions  for  a  steady  state  rather  than  at  supersaturated 
conditions. 

The  diffusion  mass  flux  ( J)  of  species  n  in  direction  §  is 


(4.1) 


where  §  is  a  dummy  variable  for  direction  x,  y  or  z. 

Energy  equation: 

V(puh)  =  V  -(kVT)-hSh.  (5) 


The  source  term  Sh  contains  contributions  from  energy  losses 
and  heat  associated  with  phase  transformations.  The  heat  source 
from  the  electrochemical  reaction  is  given  by  the  difference 
between  the  total  energy  released  by  the  electrochemical  reac¬ 
tion  at  the  cathode  membrane  surface  and  the  electrical  energy 
extracted  out  of  the  fuel  cell  [38]: 


Shr  =h 


rxn 


IA 


cv 


IF 


7Vceii  A 


CV' 


(5.1) 


The  heat  source  due  to  a  phase  change  can  be  expressed  as 


5hp  —  SW1P  x  ^fg»  (5.2) 

where  hfg  is  the  enthalpy  of  formation  of  water.  The  local  current 
density  of  the  cell  is  calculated  as  follows  from  the  open  circuit 
voltage  (Voc)  and  the  losses 

/=  —  {Voc-  Vcen-y/},  (6) 

fin 

which  can  be  rewritten  as 

Ecell  =  E0C  -TI-—I  (7) 

Cfin 

The  last  term  on  the  right  hand  side  of  Eq.  (7)  corresponds  to 
the  Ohmic  loss,  where  tm  is  the  membrane  thickness  and  crm 
is  the  membrane  conductivity,  both  of  which  are  calculated  as 
functions  of  the  water  content  on  the  membrane  surface  at  the 
anode  interface.  The  membrane  conductivity  can  be  defined  as 

0.514Mm,dryCwa  -  0.326)  exp  ( 1268  (F  _  I 
Pm, dry  )  V  \Tq  T 

(8) 


where  Tfi  =  303  K.  The  water  vapor  concentration  is  defined  as 


Pm,  dry 
37  m  ?  dry 


A. 


where  pm,dry  and  Mm? dry  are  the  material  density  and  the  equiv¬ 
alent  weight  of  a  dry  PEM,  respectively.  Water  content  in  the 
membrane  (A)  is  defined  as 


A  =  0.043  +  17. 8ax  —  39.8 a\  +  36.0 if  ax  <  1 
A  =  14  +  IA(<2k  —  1),  if  ax  >  1  ^ 


The  subscript  K  stands  for  anode  or  cathode.  aa  is  the  water 
activity  and  defined  as 


<3a  = 


3fw  ,SiP 

psat 

1  w,a 


(ID 


where  P  is  the  cell  pressure  and  Xw>a  is  the  mole  fraction  of  water 
on  anode  side. 

The  local  overpotential  (77)  for  a  PEMFC  is  included  for  both 
the  anode  overpotential  (aop)  and  cathode  overpotential  (cop), 
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and  can  be  written  as 


RT 

-  In 

ac  F 


IP 

I0q2  P0q2 


RT 

aaF 


In 


IP 

^0h2  ^0h2 


(12) 


where  P  is  the  pressure  and  Po  is  the  partial  pressure  of  the  reac¬ 
tants,  c^a  and  ac  are  the  transfer  coefficients  for  the  anode  and 
cathode,  respectively,  and  Iq  is  the  exchange  current  density.  The 
partial  pressure  is  calculated  as Po2  =  Xo2P ,  and  Pu2  =  Xu 2P. 
Here  Xq2  and  Xu2  are  the  oxygen  and  hydrogen  mole  fraction, 
respectively.  By  including  the  partial  pressure,  the  overpotential 
terms  also  include  the  activation  and  concentration  overpoten¬ 
tials. 

The  mole  fraction  of  each  species  used  in  these  model  equa¬ 
tions  is  extrapolated  to  the  membrane  surface  by  linear  extrap¬ 
olation,  except  when  a  liquid  water  film  is  generated  on  the 
membrane  surface.  In  this  latter  case,  Henry’s  Law  was  used  to 
calculate  the  solubility  of  reactants  in  the  liquid  water  film  [38]: 


~4FM° 2  =  p°2D023 


~2FMRi  =  PH2jDh2’1 


-l 


^02,glif^glif#o2,l  ^°2 


-1 


^H2,glif^glif#H2,l  *H2 


K 


a 


(13) 


where  Pgn f  is  the  pressure  at  the  gas-liquid  water  interface,  H 
the  Henry  constant  of  the  reactants  in  the  liquid  water  film  and 
tf  K  is  the  liquid  water  film  thickness  for  the  anode  and  cathode, 
which  is  defined  as 


raw,l(Emassft) 

epwiAreacv 


(14) 


where  s  is  the  porosity  of  the  GDL  and  massn  is  the  mass  of 
species  n. 

The  above  all  governing  equations  and  appropriate  boundary 
conditions  were  solved  by  using  the  user  coding  capabilities  of 
STAR-CD  that  employ  a  finite  volume  method. 


3.  Results  and  discussion 


A  series  of  simulations  were  carried  out  on  the  base  case  from 
low  operating  current  density  to  high  operating  current  density; 
the  geometrical  and  physical  parameters  for  the  base  case  are 
shown  in  Tables  1  and  2,  respectively.  The  polarization  curve  for 
the  base  case,  shown  in  Fig.  2a,  shows  three  distinct  regions  of 
losses.  The  variations  in  the  different  overpotentials  as  a  function 
of  operating  current  density,  shown  in  Fig.  2b,  also  follow  the 


Table  2 


Physical  and  electrochemical  parameters 


Parameter 

Value 

Anode  pressure  (atm) 

1 

Cathode  pressure  (atm) 

1 

Stoichiometric  rate  at  anode 

1.2 

Stoichiometric  rate  at  cathode 

2.0 

Cell  temperature  (°C) 

70 

Anode  inlet  temperature  (°C) 

80 

Cathode  inlet  temperature  (°C) 

70 

Open  circuit  voltage  (V) 

0.96 

Relative  humidity  at  anode  (%) 

100 

Relative  humidity  at  cathode  (%) 

100 

Oxygen  inlet  mole  fraction 

0.143 

Oxygen  exchange  current  density  (Am-2) 

200 

Hydrogen  exchange  current  density  (Am-2) 

2000 

Anode  transfer  coefficient 

1.2 

Cathode  transfer  coefficient 

0.6 

Porosity 

0.7 

Permeability  (m2) 

1  X  10“12 

Anode  top  surface  temperature  (z-direction)  (°C) 

70 

Cathode  top  surface  temperature  (z-direction)  (°C) 

70 

Fig.  2.  (a)  Polarization  curve  and  (b)  different  types  of  voltage  loss. 


typical  trend  found  in  the  literature  [29,39,40] .  It  should  be  noted 
that  the  Ohmic  loss  in  Fig.  2b  shows  a  non-linear  trend  rather 
than  the  linear  trend  reported  in  literature  because  the  present 
model  was  implemented  with  variable  membrane  conductivity. 


Table  1 

Geometrical  parameters 


Parameter 

Value 

Channel  length  (mm) 

34.7 

Channel  width  (mm) 

1 

Channel  height  (mm) 

1 

Membrane  length  (mm) 

31.7 

Membrane  thickness  (mm) 

0.05 

Anode  gas  diffusion  layer  (mm) 

0.25 

Cathode  gas  diffusion  layer  (mm) 

0.25 

3.1.  Effect  of  channel  cross-section  and  shoulder  width 

In  this  study,  PEMFC  configurations  with  three  channel 
cross-sections  (rectangular,  trapezoidal  and  parallelogram)  and 
different  shoulder  widths  were  used.  Cross-sectional  views  of 
these  channel  geometries  are  shown  in  Fig.  3,  and  the  geomet¬ 
rical  parameters  of  the  configurations  used  are  listed  in  Table  3. 
The  same  channel  cross-sectional  area  was  maintained,  which 
ensured  the  same  reactant  flow  rate  and  inlet  boundary  con- 


332 


D.H.  Ahmed,  H.J.  Sung  /  Journal  of  Power  Sources  162  (2006)  327-339 


b/2  a  b/2 
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Case  1 


Fig.  3.  Cross-sectional  view  of  different  geometrical  configurations. 


ditions  for  the  three  cases.  To  accentuate  any  differences  in 
cell  performance  among  the  different  channel  cross-sections, 
a  high  average  current  density  was  chosen.  The  output  voltages 
at  the  high  operating  current  density  of  /avg  =  2.4  A  cm-2  are 
shown  in  Fig.  4a  for  the  three  cases.  The  cell  voltage  for  Case 
1  (the  base  case)  is  slightly  higher  than  those  for  the  other  two 
cases.  Case  2,  which  has  a  smaller  shoulder  width  (half  shoulder 
width  =  0.35  mm)  exhibits  a  lower  cell  voltage  at  this  high  oper¬ 
ating  current  density.  To  ascertain  the  real  scenario  of  the  cell 

Table  3 


Geometrical  configurations 


Case 

Channel 
width  (mm) 

Half  shoulder 
width  (mm) 

Channel 
height  (mm) 

Reacting 
area  (cm2) 

Case  1 

0.8 

0.4 

1 

0.5076 

Case  2 

0.9 

0.35 

1 

0.5076 

Case  3 

0.8 

0.4 

1 

0.5076 

performance,  different  cell  voltage  losses  (width  average)  at  the 
membrane-GDL  interface  were  calculated  at  the  inlet  and  exit 
regions  of  the  channel  (Fig.  4b  and  c,  respectively).  Accord¬ 
ing  to  Eq.  (7),  the  cell  voltage  is  obtained  by  subtracting  the 
losses  from  the  open  circuit  voltage.  The  losses  consist  of  the 
aop,  cop  [aop  and  cop  were  already  defined  above]  and  Ohmic 
loss.  Fig.  4b  and  c  shows  that  for  all  cases,  the  aop  maintains  a 
constant  value  of  approximately  0.08  V  along  the  entire  channel 
length,  consistent  with  previous  findings  that  the  aop  is  negligi¬ 
ble  in  hydrogen  fuel  cells  [3].  By  contrast,  the  cop  and  Ohmic 
loss  account  for  around  90%  of  the  total  losses.  We  therefore 
concentrate  on  analyzing  the  Ohmic  loss  and  cop  in  detail  in  the 
following  sections. 

3.1.1.  Ohmic  loss 

According  to  Eq.  (7),  the  Ohmic  loss  is  directly  proportional 
to  the  membrane  thickness  (tm)  and  local  current  density  (/) 
and  inversely  proportional  to  the  membrane  conductivity  (am). 


Fig.  4.  (a)  Cell  voltage  obtained  for  different  cases  (/avg  =  2.4  A  cm  2),  (b)  voltage  losses  at  L/Lq  =  0.043  (entry  region)  and  (c)  voltage  losses  at  L/Lq  =  0.957  (exit 
region). 
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W/Wo 


Fig.  5.  Comparison  of  Ohmic  losses  at  the  membrane-cathode  GDL  interface, 
(a)  L/Lq  =  0.043  (entry  region)  and  (b)  L/Lq  =  0.957  (exit  region). 


Fig.  7.  Comparison  of  local  current  density  at  the  membrane-cathode  GDL 
interface,  (a)  L/Lq  =  0.043  (entry  region)  and  (b)  L/Lq  =  0.957  (exit  region). 


The  simulation  results  indicate  that  over  the  reacting  area,  the 
Ohmic  loss  is  higher  in  the  channel  area  than  in  the  shoulder 
area,  and  that  the  Ohmic  loss  decreases  on  going  from  the  chan¬ 
nel  entry  to  exit  (comparing  Fig.  5a  and  b).  The  reason  for  these 
variations  in  Ohmic  loss  in  the  PEMFC  can  be  found  in  the 
distributions  of  membrane  conductivity  (Fig.  6)  and  local  cur¬ 
rent  density  (Fig.  7).  The  membrane  conductivity  is  the  main 
factor  regulating  the  Ohmic  loss  in  the  PEMFC.  It  is  higher 
in  the  shoulder  area  and  lower  in  the  channel  area;  however, 
it  gradually  increases  along  the  channel  because  of  variations 
in  the  anode  water  activity  and  temperature.  The  distributions 
of  anode  water  activity  and  temperature  in  the  inlet  and  exit 
regions  are  shown  in  Figs.  8  and  9,  respectively.  The  temperature 
decreases  both  along  the  channel  and  away  from  the  channel  cen¬ 


ter.  The  temperature  variation  along  the  width  is  due  to  the  fixed 
temperature  and  good  thermal  conductivity  of  both  anode  and 
cathode  bipolar  plates.  The  slight  temperature  decrease  along  the 
channel  is  due  to  the  effects  of  heat  generation  from  the  elec¬ 
trochemical  reactions,  water  evaporation  and  fixed  top  surface 
temperature  of  both  anode  and  cathode  bipolar  plates.  Usu¬ 
ally,  the  membrane  water  contents  decrease  at  the  beginning 
of  the  channel  and  result  to  increase  the  temperature  sharply. 
Further  down  of  the  channel,  temperature  decreases  slowly  as 
the  membrane  water  contents  increase  through  the  electrochem¬ 
ical  reactions  [41].  These  temperature  variations  and  distribu¬ 
tions  in  Fig.  9  are  also  consistent  with  Shimpalee  and  Dutta 
[42]. 

The  above  comparison  shows  that  the  larger  channel  width 
and  smaller  shoulder  width  of  Case  2  (see  Table  3)  result  in  a 


Fig.  6.  Comparison  of  membrane  conductivity  at  the  membrane-cathode  GDL 
interface,  (a)  L/Lq  =  0.043  (entry  region)  and  (b)  L/Lq  =  0.957  (exit  region). 


Fig.  8.  Comparison  of  anode  water  activity  at  the  membrane-anode  GDL  inter¬ 
face.  (a)  L/Lq  =  0.043  (entry  region)  and  (b)  L/Lq  =  0.957  (exit  region). 
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Fig.  9.  Comparison  of  temperature  at  the  membrane-cathode  GDL  interface, 
(a)  L/Lo  =  0.043  (entry  region)  and  (b)  L/Lq  =  0.957  (exit  region). 

much  more  uniform  reaction  over  the  reacting  area.  It  is  also 
evident  that  the  temperature  at  the  membrane-cathode  GDL 
interface  is  higher  for  Case  2  than  for  Cases  1  and  3,  leading 
to  a  reduction  in  the  anode  water  activity  in  the  cell  because  the 
water  saturation  pressure  increases  with  temperature.  However, 
comparison  of  Cases  1  and  3,  which  have  the  same  shoulder  and 
channel  widths,  indicates  that  the  inclination  of  the  channel  walls 
(Case  3)  does  not  have  a  significant  effect  on  the  water  activ¬ 
ity  and  temperature  distributions,  at  this  high  operating  current 
density. 

The  results  for  the  Ohmic  loss  at  the  membrane-cathode  GDL 
interface  for  the  three  cases  (Fig.  5)  indicate  that  the  configu¬ 
ration  with  a  larger  channel  width  and  smaller  shoulder  width 
(Case  2)  experienced  higher  Ohmic  losses.  The  other  two  cases 
(Cases  1  and  3),  which  had  the  same  channel  and  shoulder 
widths,  showed  almost  the  same  Ohmic  losses. 

3.1.2.  Cathode  overpotential 

Our  results  show  that  major  losses  of  the  cell  voltage  occur  at 
the  cathode  side,  with  the  cop  accounting  for  around  65%  of  the 
total  losses  at  the  exit  region  of  the  channels.  Inspection  of  the 
cop  distribution  in  the  fuel  cell,  shown  in  Fig.  10,  discloses  that 
the  cop  increases  significantly  in  the  shoulder  region  as  well  as 
in  the  downstream  region  of  the  channels.  The  loss  due  to  the 
cop  incorporates  the  activation  and  concentration  losses  at  the 
cathode  side,  which  are  mainly  regulated  by  oxygen  availabil¬ 
ity  at  the  membrane-cathode  GDL  interface.  The  profile  of  the 
oxygen  mole  fraction  at  the  membrane-cathode  GDL  interface 
is  shown  in  Fig.  11.  Oxygen  scarcity  in  the  shoulder  region  over 
the  reacting  area  leads  to  higher  concentration  losses,  which 
become  worse  in  the  downstream  region  of  the  channel  due  to 
depletion  of  the  reactant  with  moving  downstream.  Usually,  at 
high  operating  current  density,  more  molecules  react  and  hence 
more  water  is  produced.  This  results  in  larger  quantities  of  water 
at  the  GDL,  which  reduces  the  oxygen  diffusivity  in  the  cell.  This 


effect  is  particularly  severe  in  the  shoulder  area,  where  water 
tends  to  accumulate. 

Among  the  three  cases  considered  here,  Case  2  shows  a  lower 
cop  in  both  the  channel  and  shoulder  areas.  Increasing  the  chan¬ 
nel  width  (or  face  to  the  reacting  area),  as  in  Case  2,  facilitates 
the  diffusion  of  the  reactants  to  the  reacting  area.  For  this  case, 
the  cop  is  much  lower  at  the  shoulder  region  than  for  the  other 
two  cases  (Figs.  10  and  11).  However,  Cases  1  and  3  exhibit  sim¬ 
ilar  oxygen  mole  fractions  and  cathode  overpotentials  because 
they  have  the  same  shoulder  width. 

In  the  cell  with  the  trapezoidal  cross-section  (Case  2),  the 
reactants  diffused  more  uniformly,  which  helped  to  reduce  the 
concentration  losses  and  to  obtain  a  more  uniform  local  current 
distribution  over  the  reacting  area  (see  Fig.  7).  In  the  trape¬ 
zoidal  channel  configuration  of  Case  2,  however,  the  shoulder 


Fig.  10.  Comparison  of  cathode  overpotential  at  the  membrane-cathode  GDL 
interface,  (a)  L/Lq  =  0.043  (entry  region)  and  (b)  L/Lq  =  0.957  (exit  region). 


Fig.  11.  Comparison  of  O2  mole  fraction  at  the  membrane-cathode  GDL  inter¬ 
face.  (a)  L/Lq  =  0.043  (entry  region)  and  (b)  L/Lq  =  0.957  (exit  region). 
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Table  4 


Geometrical  configurations  for  the  shoulder  width  case  studies 


Case 

Sub  case 

Ch:Sh 

Half  shoulder 
width  (mm) 

Channel 
width  (mm) 

Channel 
height  (mm) 

Reacting 
area  (cm2) 

Case  4 

Case  4.1 

1:1 

0.4 

0.8 

1 

0.5076 

Case  4.2 

1.33:1 

0.3 

0.8 

1 

0.4442 

Case  4.3 

2:1 

0.2 

0.8 

1 

0.3807 

Case  4.4 

4:1 

0.1 

0.8 

1 

0.3173 

Case  5a 

Case  5.1 

1:1 

0.4 

0.8 

1 

0.5076 

Case  5.2 

1.33:1 

0.3 

0.8 

1 

0.4442 

Case  5.3 

2:1 

0.2 

0.8 

1 

0.3807 

Case  5.4 

4:1 

0.1 

0.8 

1 

0.3173 

a  Boundary  conditions  (velocity  and  mass  flow  rate  at  anode  and  cathode  inlet)  for  Case  5  were  maintained  constant  at  the  same  values  as  for  Case  4.1. 


width  was  reduced  to  maintain  the  same  channel  cross-sectional 
area  as  for  Cases  1  and  3,  which  led  to  significantly  increased 
Ohmic  losses,  especially  in  the  shoulder  area.  The  results  of  the 
present  study  of  three  channel  geometries  indicate  that  the  shoul¬ 
der  width  is  one  of  the  major  geometrical  parameters  regulating 
voltage  losses.  Hence,  we  continued  our  study  by  examining  the 
behavior  of  cells  with  a  rectangular  channel  cross-section,  as  in 
Case  1,  but  with  different  shoulder  widths. 

In  our  investigation  of  systems  with  rectangular  channel 
cross-sections,  we  simulated  our  model  in  two  different  ways. 
The  geometrical  parameters  for  these  studies  are  shown  in 
Table  4.  For  Cases  4. 1-4.4,  the  inlet  boundary  conditions  were 
set  as  in  Eq.  (1),  and  the  shoulder  width  was  varied.  For  Cases 
5. 1-5.4,  however,  all  of  the  models  were  simulated  using  the 
same  boundary  conditions  as  for  Case  1  (i.e.,  reactant  flow  rates 
and  inlet  velocities)  but  with  different  shoulder  widths.  As  a 
result,  the  Case  5  configurations  with  smaller  shoulder  widths 
have  higher  reactant  flow  rates  and  inlet  velocities  than  for 
Case  4. 


3.2.  Channel-shoulder  width  ratio 

To  study  the  effect  of  the  channel-shoulder  width  ratio  on 
PEMFC  performance,  we  carried  out  simulations  for  different 
range  of  current  densities  like  low,  moderate  high  and  high  cur¬ 
rent  densities.  It  should  be  noted  that  in  these  simulations  we  kept 
the  channel  width  constant  and  reduced  the  shoulder  width.  Here, 
Case  4.1  is  the  same  as  Case  1  and  Cases  4.2,  4.3  and  4.4  have 
channel-shoulder  width  ratios  of  1 .33, 2  and  4,  respectively.  The 
cell  voltages  obtained  for  current  densities  of  /avg  =  2.4,  1.6  and 
0.6  A  cm-2  are  shown  in  Fig.  12a.  The  cell  voltage  for  the  oper¬ 
ating  current  density  of  0.6  A  cm-2  is  almost  independent  of  the 
channel-shoulder  width  ratio.  For  the  other  two  operating  cur¬ 
rent  densities,  however,  the  cell  voltage  varies  with  increasing 
channel-shoulder  width  ratio.  Fig.  12b  and  c  shows  the  cell  volt¬ 
age  as  a  function  of  channel-shoulder  width  ratio  for  /avg  =  2.4 
and  1.6  A  cm-2,  respectively,  with  a  smaller  voltage  scale.  It 
is  evident  in  Fig.  12b  and  c  that  the  cell  voltage  goes  through 
a  maximum  as  the  channel-shoulder  width  ratio  is  increased, 


Fig.  12.  Comparison  of  cell  voltages  for  different  channel-shoulder  width  ratios,  (a)  At  different  operating  current  densities,  (b)  /avg=2.4Acm  2  and  (c) 
/aVg  =  1.6  A  cm-2. 


336 


D.H.  Ahmed,  H.J.  Sung  /  Journal  of  Power  Sources  162  (2006)  327-339 


Fig.  13.  Comparison  of  local  current  density  distribution  for  different 
channel-shoulder  width  ratios,  (a)  Present  study  for  /avg  =  0.6  A  cm-2  at 
L/Lo  =  0.5  and  (b)  Sun  et  al.  [16]. 

indicating  that  there  exists  an  optimum  channel-shoulder  width 
ratio  of  Ch/Sh  =  1 .3-1 .4  at  high  operating  current  density.  How¬ 
ever,  the  cell  performance  decreases  significantly  for  higher 
operating  current  density  with  higher  channel-shoulder  width 
ratio.  Fig.  13a  shows  the  local  current  density  distribution 
(half  channel  and  half  shoulder)  for  different  channel-shoulder 
width  ratios  for  the  operating  current  density  of  0.6  A  cm-2  at 
L/Lo  =  0.5.  Similar  trends  of  the  local  current  density  profile 
were  observed  by  Sun  et  al.  [16]  for  operation  at  an  average 
current  density  of  around  0.65  A  cm-2.  For  Ch/Sh  =  1,  the  max¬ 
imum  local  current  density  at  the  channel-cathode  GDL  inter¬ 
face  is  just  above  the  channel-shoulder  interface  (Fig.  13a  and 
b).  On  increasing  the  channel-shoulder  width  ratio,  however, 
the  maximum  local  current  density  shifts  towards  the  shoulder 
region. 

The  different  voltage  losses  (aop,  cop  and  Ohmic  loss)  at 
the  inlet  and  exit  regions  for  systems  with  channel-shoulder 
width  ratios  of  1,  1.33,  2  and  4,  operating  at  high  current 
density  (/avg  =  2.4  A  cm-2),  are  shown  in  Fig.  14.  It  is  clear 
that  the  cop  and  Ohmic  loss  gradually  increase  and  decrease 
along  the  channel,  respectively,  for  all  of  the  channel-shoulder 
width  ratios.  At  both  the  inlet  and  exit  regions,  the  Ohmic 
loss  increases  significantly  with  increasing  channel-shoulder 
ratio.  By  contrast,  the  cop  gradually  decreases  with  increasing 
channel-shoulder  width  ratio  up  to  Ch/Sh  =  2,  and  then  remains 
almost  unchanged  on  further  increase  to  Ch/Sh  =  4.  Examination 
of  the  local  current  density  distributions  for  the  systems  with 
Ch/Sh  =  1,  1.33,  2  and  4  at  the  inlet  and  exit  regions  (Fig.  15a 
and  b,  respectively)  reveals  some  interesting  behavior.  In  the  exit 
region  (Fig.  15b),  the  local  current  density  distribution  over  the 
membrane-cathode  GDL  interface  is  more  uniform  for  higher 
channel-shoulder  width  ratios.  In  the  shoulder  area,  the  local 
current  density  is  also  improved  at  higher  channel-shoulder 
width  ratios,  however,  the  local  current  density  at  Ch/Sh  =  4  is 
not  increased  at  the  shoulder  area  compared  with  Ch/Sh  =  2. 


Fig.  14.  Comparison  of  different  voltage  losses  for  /avg  =  2.4  A  cm  2.  (a)  Entry 
region  and  (b)  exit  region. 


This  latter  observation  can  be  attributed  to  the  local  temperature 
distribution  at  the  membrane-cathode  GDL  interface.  Specif¬ 
ically,  as  shown  in  Fig.  16,  the  local  temperature  increases 
significantly  for  Ch/Sh  =  4  at  /avg  =  2.4  A  cm-2,  and  this  higher 
temperature  will  lead  to  increases  in  the  aop  and  cop  of  the  cell 
(see  Fig.  14b). 


.i  i  • 

i  !  |  i;i 

J _ 1—1 _ I _ I _ I _ 1 _ 1  1  1 _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ I _ 1 _ I _ 1 _ 1 _ I _ I _ 1 _ I _ 1  1  J _ 1 _ I _ I _ l_l _ I L 
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Fig.  15.  Comparison  of  local  current  density  distribution  for  different 
channel-shoulder  width  ratios  for  7avg  =  2.4  Acm-2.  (a)  Inlet  region  and  (b) 
exit  region. 
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Fig.  16.  Comparison  of  temperature  distribution  at  the  exit  region  for  different 
channel-shoulder  width  ratios,  (a)  7avg  =  0.6  A  cm-2  and  (b)  7av g  =  2.4  A  cm-2. 

Fig.  17  shows  the  oxygen  mole  fraction  and  membrane  con¬ 
ductivity  for  different  channel-shoulder  width  ratios  at  the  chan¬ 
nel  exit  region.  Reducing  the  shoulder  width  usually  improves 
the  oxygen  availability  at  the  reacting  area  because  the  reactants 
need  to  travel  a  shorter  distance.  This  provides  more  uniform 
reaction  in  the  fuel  cell  and,  as  a  consequence,  the  membrane 
conductivity  decreases  significantly  as  the  channel-shoulder 
width  ratio  is  increased.  A  higher  oxygen  mole  fraction  at 
the  shoulder  region  helps  to  reduce  the  cop;  however,  lower 
membrane  conductivity  increases  the  Ohmic  loss.  The  present 
simulation  results  indicate  that  the  membrane  conductivity  plays 
a  major  role  in  voltage  losses  at  high  current  density  operation 
for  higher  channel-shoulder  width  ratios. 


Membrane  width,  mm 

Fig.  17.  Comparison  of  (a)  oxygen  mole  fraction  and  (b)  membrane  conductivity 
distribution  at  the  exit  region  for  different  channel-shoulder  width  ratios  at 
7avg  =  2.4  A  cm-2. 


3.3.  Channel-shoulder  width  ratio  with  the  same  boundary 
conditions 

Above  we  examined  the  cell  performance  as  a  function  of 
channel-shoulder  width  ratio  for  systems  in  which  the  reactant 
inlet  velocities  were  defined  by  Eq.  (1).  We  now  examine  the  per¬ 
formances  of  cells  with  different  channel-shoulder  width  ratios 
(Case  5)  but  with  the  same  boundary  conditions  (reactant  mass 
flow  rates  and  inlet  boundary  conditions).  The  boundary  condi¬ 
tions  for  the  base  case  (Case  1),  which  is  the  same  as  Cases  4.1 
and  5.1,  are  used  for  the  other  sub-cases  (Case  5.2,  5.3  and  5.4). 
When  the  boundary  conditions  are  kept  constant  in  this  way, 
cells  with  a  higher  channel-shoulder  width  ratio  (e.g.,  Case  5.4) 
will  have  larger  quantities  of  reactants  and  inlet  velocities  at 
a  given  operating  current  density.  Fig.  18  shows  the  cell  volt¬ 
age  as  a  function  of  the  channel-shoulder  width  ratio  for  Cases 
4  and  5  at  an  operating  current  density  of  2.4  A  cm-2.  The  cell 
voltage  increases  significantly  at  the  optimum  channel-shoulder 
width  ratio  for  Case  5.  At  /avg  =  2.4  A  cm-2,  the  cell  voltages  for 
Cases  4  and  5  are  almost  the  same  for  channel-shoulder  width 
ratios  above  Ch/Sh  =  2.  The  different  voltage  losses  (aop,  cop 
and  Ohmic  loss)  in  the  cell  for  this  high  operating  current  den¬ 
sity  are  shown  in  Fig.  19.  Comparing  Figs.  14  and  19,  the  cop 
decreases  significantly  for  Case  5  because  higher  amounts  of 
reactants  are  flowing  through  the  channels.  However,  compared 
to  Case  4,  the  Ohmic  loss  in  Case  5  increases  to  a  much  greater 
degree  at  higher  channel-shoulder  width  ratios  because  of  the 
lower  membrane  conductivity.  Whereas  for  Case  4  there  was  no 
significant  difference  in  cop  between  Ch/Sh  =  2  and  4  (Fig.  14), 
for  Case  5  the  cop  decreases  slightly  on  going  from  Ch/Sh  =  2-4 
due  to  the  higher  amounts  of  reactants  introduced  through  the 
channels  at  the  same  inlet  temperature.  The  temperature  distri¬ 
butions  for  Ch/Sh  =  4  for  Cases  4  and  5  are  shown  in  Fig.  20.  It  is 
clear  that  the  temperature  distribution  in  the  membrane-cathode 
GDL  interface  is  lower  for  Case  5  than  for  Case  4,  which  will 
help  to  keep  a  low  cop  for  the  Case  5.  Temperature  starts  to  play 
a  critical  role  on  fuel  cell  performance  at  high  operating  cur¬ 
rent  densities  for  higher  channel-shoulder  width  ratio.  PEMFC 
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Fig.  19.  Comparison  of  different  voltage  losses  for  /avg  =  2.4  A  cm-2  by  keeping 
the  same  boundary  conditions  for  all  shoulder  widths,  (a)  Entry  region  and  (b) 
exit  region. 


with  better  cooling  systems  like  using  air  holes  (for  sub  kW  sys¬ 
tems)  or  water  holes  through  the  bipolar  plates  can  be  considered 
to  obtain  uniform  temperature  distribution  in  the  cell.  The  uni¬ 
form  temperature  distribution  will  enhance  the  cell  performance 
by  improving  the  anode  water  activity,  membrane  conductivity 
and  will  reduce  the  different  thermal  expansions  and  thermal 


stresses  of  the  membrane,  GDL  and  bipolar  materials.  However, 
the  membrane  conductivity  is  found  to  be  a  main  determinant  of 
Ohmic  loss  and  was  shown  to  be  a  major  factor  influencing  the 
total  voltage  losses  for  higher  channel-shoulder  width  ratios  at 
high  current  density  operation. 


4.  Conclusions 

In  the  present  study,  we  carried  out  simulations  of  a  PEMFC 
with  a  single  straight  channel  geometry.  The  main  aim  of  this 
work  was  to  examine  the  cell  performance  at  high  operating 
current  density  of  cells  with  different  geometric  configurations 
but  the  same  boundary  conditions.  First,  three  different  chan¬ 
nel  cross-sections  (rectangular,  trapezoidal  and  parallelogram) 
were  chosen  with  the  same  channel  height,  channel  area  and 
reacting  area,  which  ensured  the  same  boundary  conditions 
for  all  cases.  The  results  of  simulations  using  these  channel 
geometries  revealed  that  the  rectangular  channel  cross-section 
gave  higher  cell  voltages,  but  the  trapezoidal  channel  cross- 
section  gave  more  uniform  reactant  and  local  current  density 
distributions  at  the  membrane-cathode  GDF  interface.  More¬ 
over,  the  results  indicated  that  the  cathode  overpotential  and 
Ohmic  loss  were  particularly  sensitive  to  the  shoulder  width. 
We  therefore  performed  simulations  of  cells  with  a  rectangular 
channel  cross-section  and  various  shoulder  widths,  and  found 
that  the  cell  voltage  was  insensitive  to  the  shoulder  width  for 
current  densities  in  the  low  range.  At  high  operating  current 
density,  by  contrast,  there  was  an  optimum  channel-shoulder 
width  ratio  (Ch/Sh  =  1 .3-1 .4).  The  cell  performance  deteriorated 
at  higher  channel-shoulder  width  ratios,  especially  at  high  oper¬ 
ating  current  density.  Channels  with  smaller  shoulder  widths 
exhibited  higher  oxygen  concentration  and  local  current  den¬ 
sity,  but  significantly  lower  membrane  conductivity  because  of 
their  lower  anode  water  activity.  As  a  consequence,  the  cop 
decreases  and  the  Ohmic  loss  increases  in  the  fuel  cell  for  high 
operating  current  density.  The  modeling  was  repeated  so  as  to 
examine  the  variation  in  cell  performance  when  the  same  bound¬ 
ary  conditions  were  maintained  for  different  channel-shoulder 
width  ratios.  The  results  revealed  the  presence  of  an  optimum 
channel-shoulder  ratio,  at  which  the  cell  performance  increases 
significantly  for  high  operating  current  density  compared  with 
other  channel-shoulder  width  ratios.  The  membrane  conduc¬ 
tivity  and  shoulder  width  were  found  to  be  very  influential 
physical  and  geometrical  parameters,  respectively,  that  regu¬ 
lated  the  Ohmic  loss  and  overall  cell  performance  under  high 
current  density  operation. 
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